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Abstract 

We present a preliminary measurement of CP-violation parameters in the decay B° — ► K + K~ K\ 
using approximately 465 million BB events collected by the BABAR detector at SLAC. Recon- 
structing the neutral kaon as Kg — > ir + ir~ or Kg — > 7r°7r°, we analyze the Dalitz plot distribution 
and measure fractions to intermediate states. We extract CP parameters from the asymmetries 
in amplitudes and phases between B° and B° decays across the Dalitz plot. From a fit to the 
whole Dalitz plot, we measure /3 eff = 0.44 ± 0.07 ± 0.02, A C p = 0.03 ± 0.07 ± 0.02, where the 
first uncertainties are statistical and the second ones are systematic. For decays to (f>Kg, we 
measure /% = 0.13 ± 0.13 ± 0.02, A C p = 0.14 ± 0.19 ± 0.02. For decays to f K°, we measure 
/3 eff = 0.15 ± 0.13 ± 0.03, A C p = 0.01 ± 0.26 ± 0.07. From a fit to the region of the Dalitz plot with 
m K+K - > 1.1 GeV/c 2 , we measure /% = 0.52 ± 0.08 ± 0.03, A C p = 0.05 ± 0.09 ± 0.04. 
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1 INTRODUCTION 



We present a time-dependent analysis of the Dalitz plot (DP) in flavor tagged B° — > K + K~ 
decays, with the K® reconstructed as K® — ► 7r + 7r~ or — ► 7r°7r° (unless otherwise stated, charge 
conjugates are implied throughout this paper). In the Standard Model (SM), these decays are 
dominated by b — > sss gluonic penguin amplitudes, with a single weak phase. Contributions from 
b — > ugg 1 tree amplitudes, proportional to the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
V u b with a CP- violating weak phase 7 [1], are small, but may depend on the position in the Dalitz 
plot. In B° — > 4>(K + K~)K° decays the modification of the CP asymmetry due to the presence of 
suppressed tree amplitudes is at O(0. 01) [2, 3], while at higher K + K~ masses a larger contribution 
at O(O.l) is possible [4]. Therefore, to very good precision, we also expect the direct CP asymmetry 
for these decays to be small in the SM. The CP asymmetry in B° — > K + K~K® decay arises from 
the interference of decays and B° <-> B° mixing, with a relative phase of 2/3. The Unitarity Triangle 
angle (3 has been measured in B° — > [cc]K° decays to be sin 2/3 = 0.685±0.032 [5, 6]. Current direct 
measurements favor the solution of (3 = 0.37 over = 1.20 at the 98.3% C.L. [7, 8, 9, 10, 11, 12]. 
Furthermore, the (5 = 0.37 solution is the only one consistent with all indirect constraints [13, 14]. 

The decay B° — > K + K~ is one of the most promising processes with which to search for 
physics beyond the SM. Since the leading amplitudes enter only at the one-loop level, additional 
contributions from heavy non-SM particles may be of comparable size. If the amplitude from heavy 
particles has a CP-violating phase, the measured CP-violation parameters may differ from those 
expected in the SM. 

Previous BABAR measurements of the CP asymmetry in B° — ► K + K~K° decays have been 
performed on 383 x 10 6 BB events [15]. This analysis updates that previous result with a larger 
dataset. 

2 DATASET AND DETECTOR 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy B factory at SLAC. A total of 465 million BB pairs were used. 

The BABAR detector is described in detail elsewhere [16]. Charged particle (track) momenta 
are measured with a 5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) coaxial with a 1.5-T superconducting solenoidal magnet. Neutral cluster (photon) positions 
and energies are measured with an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals. Charged hadrons are identified with a detector of internally reflected Cherenkov light 
(DIRC) and specific ionization measurements (dP/dx) in the tracking detectors (DCH, SVT). 
Neutral hadrons that do not interact in the EMC are identified with detectors, up to 15 layers 
deep, in the flux return steel (IFR). 

In addition to the data collected by BABAR, this analysis uses various samples of Monte Carlo 
(MC) events based on GEANT4 [17]. A sample of simulated K + K~Kg events using a full Dalitz 
plot model based on BABAR' s previous measurement is used to study signal events, while back- 
grounds from B meson decays are studied using a separate sample of simulated events. 

3 EVENT RECONSTRUCTION 

We reconstruct B° — > K + K~Kg decays by combining two oppositely charged tracks with a K® — > 
7r + 7r~ or Kg — > 7r°7r° candidate. The K + and K~ tracks must have at least 12 measured DCH 
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coordinates, a minimum transverse momentum of 0.1 GeV/c, and must originate from the nominal 
beam spot. Tracks are identified as kaons using a likelihood ratio that combines dE/dx measured 
in the SVT and DCH with the Cherenkov angle and number of photons measured in the DIRC. 
The candidates are required to be loosely compatible with the kaon hypothesis if the K + K~ 
invariant mass is less than 1.1 GeV/c 2 , while a tighter compatibility is required in all other cases to 
further suppress background. 

For all modes, the main source of background is random combinations of particles produced in 
events of the type e + e~ — > qq (q = u,d,s,c) (continuum). Additional background from decays of 
B mesons to other final states (BB background), with and without charm particles, is estimated 
from MC events. 

We use event-shape variables, computed in the center-of-mass (CM) frame, to separate con- 
tinuum events with a jet-like topology from the more isotropic B decays. Continuum events are 
suppressed by requiring the quantity | cos 9t\ to be less than 0.9, where 6t is the angle between the 
thrust axis calculated with the B candidate's daughters and the thrust axis formed from the other 
charged and neutral particles in the event. Further discrimination comes from a Fisher discriminant 
{J-) based on 1) cos#i>) 2) 0th and 2nd order Legendre moments £i=o,2 = J2j Pj\ cos {Qj)\ l i where j 
is all tracks and clusters not used to reconstruct the B meson, pj is their momentum, and 0j is the 
angle to the B thrust axis, and 3) the magnitude of the cosine of the angle of the B with respect 
to the collision axis | cos Ob \ ■ 

In a small fraction of events, more than one B candidate in a single event passes our selection 
criteria. In this case, a single best B candidate is selected based on the K® invariant mass and on 
the quality of the kaon tracks. 

B candidates are identified using two kinematic variables that separate signal from continuum 
background. These are the beam-energy-substituted mass toes = \J {s/2 + p { ■ pb) 2 /Ef — p 2 B , 
where yfs is the total e + e~ CM energy, (Ei, p,) is the four-momentum of the initial e + e~ system and 
Pb is the B candidate momentum, both measured in the laboratory frame, and AE = Eb — \/s/2, 
where Eb is the B candidate energy in the CM frame. 

3.1 B° K+K~K° S , K° s tt+tt- 

For decays B° — > K + K~K® and K® — > it + it~ , candidates are formed from oppositely charged 
tracks with an invariant mass within 20 MeV/c 2 of the mass [1]. The vertex is required to 
be separated from the B° vertex by at least 3a. The angle oik s between the K® momentum vector 
and the vector connecting the B° and K® vertices must satisfy cos ctK s > 0.999. Distributions of 
the kinematic variables toes and AE in data, for signal and background events calculated using 
the s Vlot event-weighting method [18], are shown in Fig. 1. 

3.2 B° K K K K° s 7r°7r° 

For decays B° — ► K + K~ and K® — ► 7r°7r°, K® candidates are formed from two 7r° — ► 77 
candidates. Each of the four photons must have E^ > 0.05 GeV and have a transverse shower shape 
loosely consistent with an electromagnetic shower. Additionally, we require each 7r° candidate 
to satisfy 0.100 < m 77 < 0.155 GeV/ c 2 . The resulting K° s — ► 7r°7r° mass is required to satisfy 
0.4776 < m^o^o < 0.5276 GeV/c 2 . A K® mass constraint is then applied for the reconstruction of 
the B° candidate. 

The kinematic variables mES and AE are formed for each candidate as in Sec. 3. Distributions of 
these variables in data, for signal and background events calculated using the s Vlot event-weighting 
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m ES (GeV/c 2 ) A E (GeV) 

Figure 1: Distributions of kinematic variables (left) w-es and (right) AE for the K + K~K® (tt + tc~) 
sample. The plots show signal, with the continuum background shown in the insets. The points 
are data events weighted with the s Vlot technique, while the curves are the PDF shapes used in 
the ML fit (Sec. 4). 



method, are shown in Fig. 2. Note that the mean of the signal AE distribution is shifted from zero 
due to energy leakage in the EMC. 

4 ANALYSIS OF THE DALITZ PLOT 

Four- momentum conservation in a three-body decay gives the relation M^+m^+m^+m 2 = mf 2 + 
m 2 3 + ^23, where mfj = (pi + Pj) 2 is the square of the invariant mass of a daughter pair. This 
constraint leaves a choice of two independent Dalitz plot variables to describe the decay dynamics 
of a spin-zero particle. In this analysis we choose the K + K~ invariant mass m^+x- and the cosine 
of the helicity angle between the K + and the K® in the K + K~ center-of-mass frame, cos9h. 

We perform an extended maximum likelihood fit to the measured time dependent Dalitz plot 
distribution. We first fit on the whole DP, then fit on the mx+K- > 1-1 GeV/c 2 range (High-mass), 
then fit on the m^+K- < 1-1 GeV/ c 2 range (Low-mass). All fits are performed on the combined 
K + K~Kg(ir + ir~) and K + K~ Kg(ir°ir°) samples simultaneously. The likelihood function C for each 
subsample is defined as 



C = exp 




(1) 



where i labels the different signal and background components, j runs over all events in the sample, 
and rii is the event yield for events of the i-th. component. The probability density function (PDF) 
V% of each component is defined as 

Vi = Vi(m ES ) ■ Vi(AE) • Vi(T) ■ V D p,i{m K+K - , cos 9 H ■, At, q tag ) 8> Ki(At,<r At ), (2) 

where qt ag is the flavor of the tagged B (1 for B° and -1 for B°), and At = t rec — tt a g is the difference 
of the proper decay times of the two £>-mesons in the Y(4S) decay. <7At is the error on At, and 
1Z is the At resolution function determined from a high statistics independent sample [5]. For the 
purpose of calculating the DP coordinates m K + K - and cosOh, we refit the B candidates applying a 
B mass constraint. This ensures that the B candidates are reconstructed within the DP boundary. 
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Figure 2: Distributions of kinematic variables (left) thes and (right) AE for the K + K~ Kg(TT°Tr°) 
sample. The plots show signal, with the continuum background shown in the insets. The points 
are data events weighted with the s Vlot technique, while the curves are the PDF shapes used in 
the ML fit (Sec. 4). 

The Fisher discriminant PDF, V{J-), is only used in the Low-mass fit (see Sec. 5). Because the 
Fisher discriminant is highly correlated with the position on the DP, we do not use the Fisher 
discriminant PDF for the fit to the whole DP or for the High-mass fit. The Fisher distributions are 
shown in Fig. 3. The PDFs for the individual fit components are described in more detail below. 

4.1 Background in the Time-Dependent Dalitz Plot 

We have two background components in our fit: continuum and BB background. For the contin- 
uum background component, we use the ARGUS function [19] for "P(mEs)> an d linear polynomial 
functions for V(AE). The At distribution is described by a double-Gaussian resolution function 
convolved with a PDF of the following form: 

V(At) = f prompt 5(At) + (1 - f P rom P t)e- m/T ^, (3) 

which allows for background decays with both zero and non-zero lifetimes. The Dalitz plot for 
the continuum background is parameterized using a two-dimensional histogram PDF in the vari- 
ables m K + K - and cos9h- The histogram is filled with candidates from the region 5.2 < tties < 
5.26 GeV/c 2 . 

We estimate the amount of BB background from Monte Carlo events. The BB background is 
almost purely combinatorial and is a few percent of the total background. In the K + K~ K®(-k + tt~) 
mode, the 7tt.es and AE PDFs for the BB backgrounds are parameterized with the same functional 
forms as the continuum backgrounds. Due to non-negligible correlation between toes an d AE for 
BB background in the K + K~ Kg(ir°7r°) mode, we construct a two-dimensional smoothed histogram 
PDF in those variables. The At distribution is described with a PDF similar to the continuum 
backgrounds, but we also allow for the possibility that the non-zero lifetime component has a time- 
dependent CP asymmetry proportional to sinAm^At or cos Am^At, where Am^ is the mixing 
frequency of the B° meson. These asymmetries are set to zero in the nominal fit, but are varied as 
a systematic uncertainty. The Dalitz plot is described using a two-dimensional histogram PDF in 
a manner similar to the continuum backgrounds. 
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Figure 3: Fisher discriminant distributions for signal and continuum background and for 
K + K~Kg(ir + ir~) sample and K + K~ Kg(TT°ir ) sample. Distributions are normalized to unit area. 

4.2 Signal Decays in the Time-Dependent Dalitz Plot 

The signal components of the PDFs for "P(mEs) an d V{AE) are parameterized using modified 
Gaussian distributions: V{x) oc exp[— (x — xo) 2 /(2a± + a±(x — ^o) 2 )]- We determine the parameters 
xo, c+, u-i a+, and a_ using MC events, and fix them in fits to data. For x < xq (x > xq), the 
parameters cj_,Q!_ (a + ,a + ) are used. 

For signal events, the time-dependence is a function of location in the DP. When the flavor of 
the tagged B q tag , and the difference of the proper decay times At, are measured, the time- and 
flavor-dependent decay rate over the Dalitz plot can be written as 



dY 



1 1 e -\ At \/-r B r / _ , 

\A\ 2 + \A\ + q tag (1 - 2w) 2/m (e~ 2 ^AA*) sin Am d At 



(2?r) 3 32M| 4t b o 



q tag (l-2w) (\A\ 2 -\A\ 2 )cosAm d At], (4) 



where qtag = 1) when the other B meson is identified as a B° (B°) using a neural network 

technique [5]. The parameter uu is the fraction of events in which the B meson is mistagged with 
the incorrect flavor, and the parameter /3 is the CKM angle (3, coming from B°-B° mixing. Ap- 
proximately 75% of the signal events have tagging information and contribute to the measurement 
of CP violation parameters. After accounting for the mistag rate, the effective tagging efficiency is 
(31.2 ± 0.3)%. Events without tagging information are assigned a mistag rate of w = 0.5, and are 
included in the fit as they contribute to the determination of the Dalitz plot parameters. Decay 
amplitudes A and A are defined in (6) and (7) below. M B o, t b o, and Am^ are the mass, lifetime, 
and mixing frequency of the B° meson, respectively [1]. 
The PDF for the Dalitz plot rate takes the form 

V DP oc dT (m K+K -, cos 6 H , At, q tag ) ■ e(m K+K -,cos6 H ) ■ |J| ®K(At,a A t), (5) 
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where \J(m,K+K-)\ = (^ m K+K-)(^QP) 1S the Jacobian of the transformation (m 2 K + K --,rn 2 K+K a) 
{m K + K -, cos 9h), and is given in terms of the charged kaon momentum q and neutral kaon momen- 
tum p, in the K + K~ frame. The efficiency e is calculated from high-statistics samples of simulated 
events and depends on the position on the Dalitz plot. 

The amplitude A (A) for the decay B° -> K + K~K° (B° -> K~K+K0) is, in our isobar model, 
written as a sum of decays through intermediate resonances: 

A = ^c r (l + & r )e i( ^ +<v) -f r , and (6) 

r 

A = ^c r (l-6 r )e^-^-/ r . (7) 

r 

The isobar coefficients c r and r are the magnitude and phase of the amplitude of component r, and 
we allow for different isobar coefficients for B° and B° decays through the asymmetry parameters 
b r and S r . The function f r = F r xT r x Z r describes the dynamic properties of a resonance r, where 
F r is the form-factor for the resonance decay vertex, T r is the resonant mass-lineshape, and Z r 
describes the angular distribution in the decay [20, 21]. 

Our model includes the 0(1020), for which we use the Blatt-Weisskopf centrifugal barrier factor 
F r = l/\/l + (Rq) 2 [20], where q is the daughter momentum in the resonance frame, and R is the 
effective meson radius, taken to be R = 1.5 GeV -1 (0.3 fm). For the scalar decays included in our 
model (/o(980), X(1550), and Xco), we use a constant form-factor. Note that we have omitted a 
similar centrifugal factor for the B° decay vertex into the cpK intermediate state since its effect is 
negligible due to the small width of the 0(1020) resonance. 

The angular distribution is constant for scalar decays, whereas for vector decays Z ~ q ■ p, 
where q is the momentum of the resonant daughter, and p is the momentum of the third particle 
in the resonance frame. We describe the line-shape for the 0(1020), X(1550), and XcO using the 
relativistic Breit-Wigner function 

T{m) = — 2 " ■ r , v ( 8 ) 

mf. — rn K + K - — im r l (m) 

where m r is the resonance pole mass. The mass-dependent width is given as T(i71k+k-) = 
T r {q/q r ) 2L+1 irnr/f^K+K-) (Fr{q)/ F r (q r )) 2 , where L is the resonance spin and q = q r when m K + K - 
m r . For the 0(1020) and XcO parameters, we use average measurements [1]. The X(1550) is less 
well-established. Previous Dalitz plot analyses of B + -> K+K+R- [22, 24] and B° -> K + K~K° 
decays [25] report observations of a scalar resonance at around 1.5 GeV/c 2 . The scalar nature has 
been confirmed by partial- wave analyses [23, 24]. However, previous measurements report inconsis- 
tent resonant widths: 0.145±0.029 GeV/c 2 [22] and 0.257± 0.033 GeV/c 2 [24]. Branching fractions 
also disagree, so the nature of this component is still unclear [26]. In our nominal fit, we take the 
resonance parameters from Ref. [24], which is based on a larger sample of BB decays than Ref. [22], 
and consider the narrower width given in the latter in the systematic error studies. 
The /o(980) resonance is described with the coupled-channel (Flatte) function 

T(m K+K -) = — 2 ——( 1 V' ( 9 ) 



where PK{m K+K -) = 2yl - 4m 2 K /m 2 K+K _, p^{m K+K -) = 2yl - 4m 2 r /m 2 K+K _, and the coupling 

strengths for the KK and -kit channels are taken as g w = 0.165±0.018 GeV/c 2 , gx/g-K = 4.21±0.33, 
and m r = 0.965 ± 0.010 GeV/c 2 [27]. 
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In addition to resonant decays, we include non-resonant amplitudes. Existing models consider 
contributions from contact terms or higher-resonance tails [28, 29, 4], but they do not capture 
features observed in data. We rely on a phenomenological parameterization [22] and describe the 
non-resonant terms as 

A NR (Anr) = (c 12 e^ 12 e- am ? 2 + c 13 e^ 13 e" Qm » + C2 ^e~ am ^) ■ (1 ± b NR ) ■ e ±lS ™, (10) 

where 1,2,3 denote the three daughter particles of the B meson. The slope of the exponential 
function is consistent among previous measurements in both neutral and charged B decays into 
three kaons [22, 24, 25], and we use a = 0.14 ± 0.02 GeV~ 2 • c 4 . 

We compute the direct CP-asymmetry parameters for resonance r from the asymmetries in 
amplitudes (b r ) and phases (8 r ) given in Eqs. (6, 7). We define the rate asymmetry as 

and (3 e ff(r) = j3 + 5 r is defined as the total phase asymmetry. These asymmetries are related to the 
CP asymmetry parameters C and —rjS using the approximations 

C r « -A C p(r), and (12) 
-VrS r ^^^sm(2p eff (r)), (13) 

where r\ r is the CP eigenvalue of the final state. The fraction for resonance r is computed as 

jdcos9 H dm K+K - ■ \J\ • (|A| 2 + IA-1 2 ) (u) 
Jdcos 6 H dm K+K - ■ \J\ ■ (\A\ 2 + \A\ 2 ) ' 1 ' 

The sum of the fractions can be different from unity due to interference between the isobars. 

In addition to the previously mentioned resonances, the decays B° — > D + K~ (D + — > K + K®) 
and B° — ► DfK~ (Df — > K + K®) are also counted as signal. We include non-interfering amplitudes 
for these modes in our Dalitz plot model, parameterizing the Dr s \ mesons on the Dalitz plot as 
Gaussian distributions with widths taken from studies of simulated events. The parameters b r and 
S r are fixed to zero for the decays B° — > D + K~, B° — > DfK~, and B° — > XcoK® throughout this 
analysis. 

5 RESULTS 

In order to determine parameters of the Dalitz plot model, we perform three fits: 1) whole DP fit, 
2) Low-mass (m K + K - < l.lGeV/c 2 ) region fit, and 3) High-mass {m K + K - > l.lGeV/c 2 ) region 
fit. 

5.1 The whole Dalitz Plot fit 

We perform a fit to both 4316 B° -► K + K~ K Q s (ir+ n~) and 2205 B° ^+K-K0(7r 7r°) candidates 
simultaneously in the full Dalitz plot. In this step we assume that all charmless decays have the 
same CP-asymmetry parameters. A Fisher discriminant cut (—2.5 < T < 4), which retains about 
95% of signal events and 60% of continuum events, is applied. We do not include the Fisher PDF 
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in the fit. We vary the event yields, isobar coefficients, and the two CP-asymmetry parameters 
Aqp and (3 e g averaged over the Dalitz plot. We find a signal yield of 1268 ±43 (K + K~ K®(ir + 
and 160 ±19 (K + R- K°(it it )) events, and a BB background yield of 47 ±31 {K + K~ K°(tt + 'tt ■ )) 
and 24 ± 16 (K + K~ Kg(ir°7r )) events. The isobar amplitudes, phases, and fractions are listed in 
Table 1. The resonant fractions do not add up to 100% due to interference between the resonances. 
The CP-asymmetry parameters, and the correlation coefficients p between them, are summarized 
in Table 2. Fig. 4 shows a projection of the Dalitz plot variable m K + K - . Fig. 5 shows distributions 
of At for i?°-tagged and _B°-tagged events, and the asymmetry A(At) = (N B o —N^ )/(N B o +N^ ), 
obtained with the s Vlot technique. 

To calculate the significance of the nominal f3 e ff result, many fits are performed with fixed but 
different (3 e g values. The change in likelihood as a function of fi e g is shown in Fig. 6. 



Decay 


Amplitude c r 


Phase 4> r 


Fraction T r (%) 


</>(1020)i^ 


0.00897 ± 0.00096 


-0.341 ±0.232 


12.6 ± 1.0 


/o(980)tf° 


0.542 ± 0.044 


-0.201 ±0.157 


27.8 ±7.1 


X (1550)K° 


0.141 ±0.017 


-0.370 ±0.154 


5.70 ± 1.70 


NR {K+K-) 


1 (fixed) 


(fixed) 


98.1 ± 18.7 




0.328 ± 0.058 


1.81 ±0.23 


10.5 ±3.4 


(K-K° s ) 


0.353 ± 0.066 


-1.44 ±0.27 


12.1 ±3.8 


XcoK° s 


0.0298 ± 0.0046 


0.732 ± 0.437 


2.53 ±0.60 


D+K- 


1.34 ±0.19 




3.43 ± 0.69 


DfK- 


0.826 ±0.160 




1.37 ±0.46 



Table 1: Isobar amplitudes and phases from the fit to the whole DP. Three rows for non-resonant 
(NR) contribution correspond to coefficients of exponential functions in Eq. (10). The errors are 
statistical only. 



Name 


Whole DP 


High-mass 


A C p 


0.03 ± 0.07 ± 0.02 


0.05 ± 0.09 ± 0.04 


/% 


0.44 ±0.07 ±0.02 


0.52 ±0.08 ±0.03 


P 


0.041 


0.031 



Table 2: The CP-asymmetry parameters from the whole DP fit and the High-mass region fit. The 
first error is statistical and the second is systematic, p is the correlation coefficient. 

5.2 High-mass fit 

We perform a fit to both 3112 B° -► K+ K~ K^+tt-) and 1917 B° -► K+ K~ K°(n°ir ) candidates 
in the High-mass region {m^+K- > l.lGeV/c 2 ) simultaneously. We fix all isobar coefficients to 
the values from the whole DP fit. We vary yields and shared CP- asymmetry parameters. We 
find a signal yield of 894 ± 36 K+ K~ K° s (it + it-) and 117 ± 16 K~ K°(tt c V ) events, and a BB 
background yield of 50 ± 31 {K+ K~ K° s {n + ir-)) and 20 ± 15 (i^+^-^(7r vr )) events. The fit 
results are summarized in Table 2. Fig. 7 shows a projection of the Dalitz plot variable cos 9 B for 
events in this region, using the s Vlot technique. 
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Figure 4: For the whole DP region fit, the distribution of the Dalitz plot variable m^+K- f° r 
signal-weighted data events (points) compared with the fit PDF in the K + K~ K®(tt + tt~) mode. 
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Figure 5: The At (top) distributions and asymmetries (bottom) in the whole DP (left) and Low- 
mass region (right), for the K + K~ K®(tt + it~) mode. For the At distributions, B°- (B -) tagged 
signal- weighted events are shown as filled (open) circles, with the PDF projection in dashed red 
(solid blue). 
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Figure 6: The change in the value of -21og(£) as a function of [3 e g, for (a) the whole DP, (b) the 
High-mass region, (c) / (980), and (d) 0(1020). 
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Figure 7: For the High-mass region fit, the distribution of the Dalitz plot variable cos Oh for 
signal-weighted data events (points) compared with the fit PDF in the K + K~ K®{ir + mode. 



5.3 Low-mass fit 

In order to measure CP-asymmetry parameters for components with \ow-K + K~ mass with reduced 
model-dependence from the rest of the DP, we select events using a cut of m K + K - < 1.1 GeV/c 2 . 
Because we are only selecting a small region of the DP, the correlation between the Fisher discrimi- 
nant T and the DP location is unimportant. We therefore relax the cut on and add the T PDF 
to the fit. After these requirements on m^+K- an d J~, there are 1846 (K + K~ K®(ir + ir~)) and 493 
(K + K~ Kg(TT°ir )) candidates remaining. The most significant contributions in this region come 
from (p(W20)Kg and fo(980)K® decays, with a smaller contribution from a low-K + K~ mass tail of 
non-resonant decays. We fix all the isobar coefficients except for those of the 0(1020) to the values 
from the whole DP fit, and fix the CP-asymmetry parameters b r and S r for all resonances except 
the 0(1020) and /o(980) to be 0. We vary the events yields, isobar coefficients for the 0(1020), and 
separate CP-asymmetry parameters for the 0(1020) and /o(980) in the fit. We find signal yields of 
381 ± 23 (K+K-K%(ir + ir-)) and 40 ± 9 (K + K~ K°(7r°vr )) events, and BB background yields of 
12 ± 13 (K + K-Kf(ir + ir-)) and -3 ± 5 (K+ K~ K°(7r°7r )) events 

The CP-asymmetry results are listed in Table 3; the systematic uncertainties will be described 
in Sec. 6. We find two solutions with likelihood difference Alog(£) = 0.1. Solution (1) is consistent 
with the SM, while Solution (2) has a value of (5 e g for the /o(980)i^^ decay that differs significantly 
from the SM, as shown in Table 3. The two solutions also have significantly different values of c r 
for the 0(1020). Both solutions also have a mathematical ambiguity of ±7r radians on [3 e ff for the 
0(1020), and a correlated ambiguity of ±7r radians on the isobar parameter r for the 0(1020). This 
ambiguity is present because the decay amplitude contains interference terms that only depend on 
the linear combinations f3 e ff + 4> T and (5 e jj — r . We choose Solution (1) as our nominal solution. The 
correlation coefficients p between the CP parameters for Solution (1) are shown in Table 3. Because 
the decay rate depends on interference terms between the 0(1020)-?^ and /o(980)i^T^ decays, the 
significant correlation between the measured CP parameters is expected. 

Fig. 8 shows projections of the Dalitz plot distributions of events in this region, using the s Vlot 
technique. Fig. 5 shows distributions of At for l?°-tagged and -B°-tagged events, and the asymmetry 
A(At) = (N B o-N Eo )/(N BO +N Eo ). 

The decay B° — > 0(1020)7^, with highly suppressed tree amplitudes, is, in terms of theoretical 
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Name 


Solution (1) 


Solution (2) 


Correlation 

1 9 A 


1 Acp^K") 
2(3 eff (^K Q s ) 
3 A CP (f K° s ) 

iPeff(f K° s ) 


0.14 ±0.19 ±0.02 
0.13 ±0.13 ±0.02 
0.01 ±0.26 ±0.07 
0.15 ±0.13 ±0.03 


0.13 ±0.18 
0.14 ±0.14 
-0.49 ± 0.25 
3.44 ±0.19 


1.0 -0.09 -0.28 0.09 
1.0 0.54 0.65 
1.0 0.25 
1.0 



Table 3: CP-violation parameters for B° — » K + K for tn K + K - < 1.1 GeV/c 2 . The first error 
is statistical and the second is systematic. Correlation coefficients are given for Solution (1) only. 
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Figure 8: For the Low-mass region fit, the distributions of the Dalitz plot variables rnx+K- (left) 
and cos 9h (right) for signal-weighted data events (points) compared with the fit PDF in the 
K + K-K^(tt + tt-) mode. 



uncertainty, the cleanest channel to interpret possible deviations of the CP-violation parameters 
from the SM expectations. Values of (5 e g are consistent with the value found in [cc]-fT° decays [5, 6]. 

We also calculate the parameters C and — tjS for ^(1020)^ and fo(980)K® using the expressions 
in (12) and (13). The results are shown in Table 4, along with C and — rjS for the whole DP and 
High-mass fits. 

6 SYSTEMATIC STUDIES 

We study systematic effects on the CP-asymmetry parameters due to fixed parameters in the uies 
and AE PDFs. We assign systematic errors by comparing the fit with nominal parameters and with 
parameters varied by their error (±1<t), and assign the average difference as the systematic error. 
In addition, we account for a potential fit bias using values observed in studies with MC samples 
generated with the nominal Dalitz plot model. We take the average values of the bias observed in 
these studies as the systematic error. We account for fixed At resolution parameters, B° lifetime, 
mixing and flavor tagging parameters. We also assign an error due to interference between 
the CKM-suppressed b — ► ucd and the favored b — > cud amplitude for some tag-side B decays [32] . 
Smaller errors due to beam-spot position uncertainty, detector alignment, and the boost correction 
are based on studies done in charmonium decays. In the cases of the Low- mass and High-mass 
fits, we also assign systematic errors due to the isobar coefficients that are fixed to the result from 
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C 




Whole DP 


-0.03 ±0.07 ±0.02 


0.77 ±0.09 ±0.02 


High-mass 


-0.05 ±0.09 ±0.04 


0.86 ±0.08 ±0.03 


0(1020)^ 
/o(980)^ s ° 


-0.14 ±0.19 ±0.02 
-0.01 ±0.26 ±0.07 


0.26 ±0.26 ±0.03 
0.29 ±0.25 ±0.06 



Table 4: The CP asymmetry parameters C and — rjS, derived using Equations (12) and (13). 
Results are shown for the whole DP, the High-mass region, and for both (p(1020)Kg and fo(980)Kg 
in the Low-mass region. For the Low-mass results, only Solution (1) is shown. The first error is 
statistical and the second is systematic. 

the whole DP fit. In all fits we assume no direct CP violation in decays dominated by the b — > c 
transition {xcqK% D^K). 

We also assign an error due to uncertainty in the resonant and non-resonant line-shape parame- 
ters. The systematic uncertainty associated with the resonant component includes the uncertainty 
in the mass and width of the X(1550), estimated by replacing the parameters used in the nominal 
fit with the values found by different measurements: m r = 1.491 GeV/c 2 , T = 0.145 GeV [22]. All 
the systematic uncertainties are summarized in Table. 5. 



Parameter 


Whole DP 

A CP Peff 


« fo^> 

AcP Peff A C p I3 e ff 


High-mass 

A CP Peff 


Fixed PDF Parameters 
Fit Bias 

DCSD, Beam Spot, other 
Dalitz Model 


0.010 0.010 
0.007 0.011 
0.015 0.004 
0.005 0.005 


0.014 0.010 0.025 0.015 
0.009 0.012 0.011 0.011 
0.015 0.004 0.015 0.004 
0.009 0.002 0.060 0.024 


0.013 0.010 
0.014 0.009 
0.015 0.004 
0.027 0.023 


Total 


0.020 0.016 


0.024 0.016 0.068 0.031 


0.036 0.026 



Table 5: Summary of systematic errors on CP-asymmetry parameters. Errors for (frK® and foK® 
CP-parameters are based on the Low-mass region fit. Total is obtained from the quadratic sum of 
the individual systematics. 



7 CONCLUSIONS 

We performed a ML fit to analyze the DP distribution of B° — ► K + K~K® decay with the full BABAR 
dataset. From a fit to the whole DP, we measure /% = 0.44 ±0.07 ±0.02, A C p = 0.03 ±0.07 ±0.02, 
consistent with our previous measurements [15] and compatible with the Standard Model values 
(3 ~ 0.37, Aqp = 0. We measure CP violation with a significance of 6.7 standard deviations 
(including statistical and systematic errors), and we reject the solution near tt/2 — /3 at 4.8 standard 
deviations. 

From a fit to the region of the DP with m K + K - > 1.1 GeV/ c 2 , we measure f3 e jj = 0.52 ± 0.08 ± 
0.03 and Acp = 0.05 ±0.09 ±0.04, compatible with the Standard Model expectations. We measure 
CP violation in this High-mass region at 6.7 standard deviations. 

From a fit to events at low K + K~ masses, we measure (3 e jj = 0.13 ± 0.13 ± 0.02 and Aqp = 
0.14 ±0.19 ±0.02 for B° -► ^(1020)^°, and /% = 0.15 ± 0.13 ± 0.03 and A CP = 0.01 ± 0.26 ± 0.07 
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for B° — ► foKg. The results for /3 ej ff are roughly 1.7 standard deviations below the Standard Model 
value. 

These results supersede our previous measurements [15] made on a smaller dataset. All of our 
results are consistent with our previous measurements. 
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